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Abstract 

Wetting phenomena plays an interesting role in the technological development of materials. Recently, 

much attention has been directed to the study of magnetic solid films. To understand, theoretically, the effect 

of surface on wetting and layering transitions in these systems, we will give in this chapter a presentation of 

several relevant recent theoretical works realized on wetting and layering transitions in which we explain how 

the nature of surface can affect the behavior of wetting and layering transitions. Using different spin systems 

models and different numerical and approximate methods such as mean field and effective field theories, real 

space renormalization group technique, transfer matrix method and Monte-Carlo simulations, it is found that 

the wetting and layering transitions depend on the nature of the surface magnetic field, the surface coupling 

strength, the surface crystal field, the geometry of surface, the in-homogeneity of substrate and the quantum 

fiuctuations at the surface. Indeed the layering transition temperature increases or decreases as a function of 

the thickness of the film depending on which the surface coupling is smaller or greater than a critical value. 

Thus the wetting transition temperature can be greater or smaller than the layering temperature. In the case 

of variable surface crystal field an order-disorder layering reentrance appears with new critical exponents 

and multi-critical behaviors. Wetting phenomena is also influenced by the curvature of the surface. In 

fact the intra-layers transitions appear under the effect of edge on the surface and/or corrugated surface. 

However, there exist an intra-layering surface temperature above which the surface and the intra-layer surface 

transitions occur . While the bulk layering and intra-layering transitions above an other finite temperature 

which is greater than the former one. Such behavior is not seen in the case of a perfect surface. Beside 

this and under the corrugation effect , the pre-wetting phenomena and layering transitions occur at zero 

temperature, even in the case of a uniform magnetic field. In this case the wetting temperature depends 

strongly on the width of the corrugated-steps but weakly on its number. The quantum fiuctuations introduce 

a new wetting transitions, namely the wetting, layering and roughening transverse fields. These fiuctuations 

lead to the appearance of a sequence of layering sublimations transitions which occur before bulk melting 

above a critical transverse field which depends on temperature and/or chemical potential. Moreover, the 

wetting temperature depends strongly on the randomness of the quantum fluctuations. The chapter is 

organized as follows: Section 1 is reserved to the Introduction and deflnitions. In section 2 we present 
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the relation between the nature of surface and wetting phenomena. The effect of the surface geometry on 
wetting and layering transitions is presented in section 3. While quantum fluctuations effects are given in 
section 4. References are cited in section 5. 
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1 Introduction 



Recently, the wetting and layering transitions of magnetic Ising systems have been studied by several authors. 
Multilayer films adsorbed on attractive substrates may exhibit a variety of possible phase transitions, as has been 
reviewed by Pandit et al. [1] , Pandit and Wortis [2] , Nightingale et al. [3] and Ebner et al. [4-7] . A variation of 
phase diagrams with the strength of the substrate potential in lattice gas model for multilayer adsorption has 
been studied by Patrykiejew et al. [8] using Monte Carlo simulations and molecular field approximation. One 
type of transitions is the layering transitions, in which the thickness of a solid film increases discontinuously by 
one layer as the pressure is increased. Such transitions have been observed in a variety of systems including for 
example ^He [9,10] and ethylene [11,12] adsorbed on graphite. 

Many experiments works have for subject the study of wetting transitions. However, the effect of substrate 
roughness on wetting has been studied [13,14] and the thickness profile of the drop of a nematic liquid crystal 
on a rough surface [15]. The effect of the nature of the adsorbate has been studied by several authors namely 
[16] the critical wetting of pentane on water, the critical adsorption and the layering transitions, at the free 
surface, of a smectic liquid crystal [17] , the effect of segregation on wetting in films of a partially miscible polymer 
blend[18] , the effect of the chain length of alkanes on water [19] , the critical adsorption and dimensional crossover 
in epitaxial FeCo films [20], the heterogeneous hole nucleation in electron-charged '^He wetting films [21], the 
existence of tricritical wetting in liquid mixtures of methanol and the n-alkanes [22] , the microscopic structure of 
the wetting film at the surface of liquid Ga-Bi alloys [23] , the tetra point wetting at the free surface of liquid Ga- 
Bi, in which a continuous surface wetting transition, pinned to a solid-liquid-liquid-vapor tetra coexistence point 
occurs [24], the short range wetting at liquid (Ga-Bi) alloy surfaces [25], the critical Casimir effect and wetting by 
Helium mixtures [26]. The dynamic aspect has been studied by several authors; namely, the dynamic properties 
of a fluid interface of nanometric curvature formed by capillary condensation between solid substrates [27] , the 
dynamic of wetting layer [28], and the dynamic of layering transition in confined liquids [29]. The nature of the 
transitions has attracted attention like; the wetting near the triple point [30] , and the sequence of two wetting 
transitions induced by tuning the Hamaker constant [31]. However, a simple lattice gas model with layering 
transitions and critical points has been introduced and studied in the mean field approximation by de Oliveira 
and Griffiths [32]. These experimental studies have motivated many theoretical works. However several methods 
have been used to study wetting and layering transitions. The transfer matrix and scaling theory have been 
used to study the curvature controlled wetting in two dimensions [33] , the adsorption of polyelectrolytes at an 
oppositely charged surface [34] . Using Monte Carlo simulations the effects of van der Waals surface interactions 
on wetting transitions in polymer blends [35], the wetting properties of neon on heterogeneous surfaces [36], 
the wetting transitions of hydrogen and deuterium on the surface of alkali [37], the multilayer adsorption of 
binary mixtures [38], the competition between short-ranged attractions and long-ranged repulsions in systems 
like langmuir monolayers, magnetic films, and adsorbed monolayers using numerical simulations and analytic 
theory [39] , the double wedge wetting [40] , and the wetting behavior of associating binary mixtures at attractive 
walls [41]. The variational and mean field approximation have been used to investigate many interest problems 
such as; the dewetting, partial wetting and spreading of a two-dimensional monolayer on solid surface [42], the 
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line tension between fluid phases and a substrate [43] , the effect of confinement and surface enhancement on 
superconductivity [44] , the symmetric polymer blend confined into a film with antisymmetric surfaces [45] , the 
wetting of potassium surfaces by superfluid ^iJe [46] , the filling transition for edge [47] , the wetting of film on 
chemically heterogeneous substrates [48], the wetting and the interaction potential between spherical particles 
[49] , the phase diagram for morphological transitions of wetting films on chemically structured substrates [50] , 
the phase behavior of colloidal binary hard-platelet and hard-rod fiuids [51], the influence of random self-afiine 
and mound substrate roughness on the wetting of adsorbed van der Waals films [52] , the wetting of van der Waals 
solid films on self-affine rough surfaces [53] , the preroughening and reentrant layering transitions on triangular 
lattice substrates [54] , the layering transitions, disordered flat phases, reconstruction, and roughening of a semi- 
infinite crystal [55]. Other method have been used to study wetting phenomena such as; the renormalization 
group of the capillary parameter at complete wetting [56], the surface deconstruction and roughening in the 
multiziggurat model of wetting [57], the morphological transitions of wetting layers on structured surfaces 
[58], the limiting height at which a solid can approach an air-water surface without becoming wet [59], the 
surface roughness and hydrodynamic boundary slip of a newtonian fluid in completely wetting system [60], 
the modeling at low temperature of a ^iJe film adsorbed on a Li substrate [61] and the prewetting of "^He 
on rubidium [62], the critical adsorption on defects in Ising magnets and binary alloys [63], the exact study 
of edge filling transition for Ising corners [64], the wetting hysteresis at the molecular scale [65], the wetting 
transition in a two-dimensional Ising model with a free rough surface [66]. Ebner and Saam [67] carried out 
Monte Carlo simulations of such a lattice gas model. Huse [68] applied renormalizaion group technique to this 
model. It allowed the study of the effects on an atomic scale or order disorder transitions in the adsorbed layers, 
which may have considerable influence on the layering transitions and tracing back macroscopic phenomena on 
inter-atomic potentials. Benyoussef and Ez-Zahraouy have studied the layering transitions of Ising model thin 
fllms using a real space renormalization group [69], and transfer matrix methods [70]. Using the mean field 
theory. Hong [71], have found that depending on the values of the exchange integrals near the surface region, 
the film critical temperature may be lower, higher than, or equal to that of the bulk. 

On the other hand, the development of molecular-beam epitaxy has allowed the fabrication of high quality 
thin films [72]. Among them, the magnetic film consisting on a few atomic layers is of particular importance, 
as it has been used as a convenient model for theoretical and experimental analysis of a variety of magnetic 
phenomena. One other problem investigated is the critical temperature at which the system undergoes a 
magnetic-nonmagnetic phase transition. Of particular interest is the change of this temperature when the 
magnetic system becomes finite in a given direction. 

Briefly, depending on the values of magnetic coupling constants near the surface region, two different behav- 
iors at the surface can be demonstrated [71] : (i) the surface critical temperature is the same as of the bulk one 
when the coupling constant on the surface is smaller than a critical value "ordinary transition" . (ii) The surface 
critical temperature may be higher than the corresponding bulk one when the surface coupling is larger than 
the critical value "extraordinary transition". At the critical ratio, one encounters the "special transition point", 
with critical properties of the surface transition deviating from those at the ordinary or the distinct surface 
transition. As for the film system, which is finite in one direction, it has been established that its magnetic 
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properties can differ greatly from those of the corresponding bulk [73-77] . Experimental results [78] showed that 
the critical temperature of a vanadium film depends on the film thickness and its critical behavior is like that 
of the two-dimensional system rather than that of the three-dimensional bulk. 

Since the theory of surface critical phenomena started developing, much attention has been devoted to 
the study of the Blume Capel (BC) model over semi-infinite lattices, with modified surface couplings. Phase 
diagrams in the mean field approximation reporting four possible topologies at fixed bulk/surface coupling ratios 
have been determined [79-83]. An analogous analysis have also been done using a real space renormalization 
group transformation [84] . Some other works referring to particular regions of the phase space are, for example, 
those using: the mean field approximation [85] , the effective field approximation [86] and the low temperature 
expansion [87] . These works show that it is possible to have a phase with ordered surface and disordered bulk, 
which is separated from the completely ordered phase by the so-called extraordinary transition and from the 
completely disordered phase by the surface transition. In the absence of this phase, the transition between the 
completely ordered and the completely disordered phase is called ordinary. These three kinds of phase transitions 
have a meeting point named special which is generally a multi-critical point. The discussion presented in Ref. 
[88] shows that the strong interest in these models arises partly from the unusually rich phase transition behavior 
they display as their interaction parameters are varied, and partly from their many possible applications. The 
bilinear interaction considered in most of these cases is ferromagnetic. The spin-1 Ising systems are used to 
describe both the order-disorder transition and the crystallization of the binary allow, and it was solved for 
the anti- ferromagnetic case in the mean field approach [89]. The decomposition of a line of tri-critical points 
into a line of critical end points and one of double critical points is one of the most interesting features of the 
mean field phase diagram for the anti- ferromagnetic spin-1 Blume-Capel model in an external magnetic field 
[90]. The transfer-matrix and Monte Carlo finite-size-scaling methods [91], are also applied to study this model 
but such decomposition does not occur in the two dimensional case. The finite cluster approximation has been 
applied by Benyoussef et al. [92] in order to study the spin-1 Ising model with a random crystal field. On the 
other hand, the transverse field or crystal field effects of spin-1 Ising model has been studied by several authors 
[93-96]. 

The experimental measurements of layer-by-layer ordering phenomena have been established on free-standing 
liquid crystals films such as nmOBC (n — alkyl — 4' — n — alkyloxybiphenyl — 4 — carboxylate) [97,98] and 
biCOOBC ( n-pentyl-4'-n-pentanoyloxy-biphenyl-4-carboxylate) [99] for several molecular layers. More re- 
cently, Lin el al. [100] have used the three-level Potts model to show the existence of layer-by-layer ordering 
of ultra thin liquid crystal films of free-standing 54:COOBC films, by adjusting the interlayer and intra layer 
couplings between nearest- neighboring molecules. 

Depending on the value of the crystal field at the surface some systems undergo a phase transition from a 
completely disordered state, at high temperatures, to a state with the first k layers disordered and the rest of 
the layers ordered at low very temperatures. One can then observe a particular behavior when the transition to 
the order state of the fc*'' layer occurs not only with the decrease but also with the increase of the temperature. 
Such phase transitions are called the reentrant phenomena. 

The reentrant first-order layering transitions have been observed experimentally in multilayer argon films 
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on graphite by Youn and Hess [101]. Interesting results were obtained within MFA such [102] as order-disorder 
layering transitions and reentrant phenomena. But it is well known that the MFA approximation does not take 
into account the thermal fluctuations and neglects correlations between spins. This can lead to incorrect results. 

The Blume-Capel model (BC) was originally proposed to study the first-order magnetic phase transitions 
in spin— 1 Ising systems [103]. This model was generalized to the Blume-Emery-Grifiiths (BEG) to study 
phase separation and superfluity in ^He-*He mixtures [104]. Later it has been applied to describe properties 
of multicomponent fluids [105], semiconductor alloys [106] and electronic conduction models [107]. The (BC) 
model is not exactly solvable in more than one dimension, but it has been studied over inflnite d-dimensional 
lattices by means of many different approximate techniques and its phase diagram is well known. 

On the other hand, ferroelectric films can be described by an Ising model and when the film becomes 
very thick, its properties are those of the semi-inflnite Ising system [108-110]. From the experimental point 
of view, the most commonly studied magnetic multilayers are those of ferromagnetic transition metal such as 
Fe/Ni, where the coupling can exist between magnetic layers [111-113]. The discovery of enormous values of 
magnetoresistance in magnetic multilayers are far exceeding those found in single layer films and so exceeds 
the discovery of oscillatory interlayer coupling in transition metal multilayers. These experimental studies have 
motivated much theoretical works to study magnetic thin films as well as critical phenomena [114-117]. This 
is partly motivated by the development of new growth and characterization techniques, but perhaps more so 
by the discovery of many exciting new properties, some quite unanticipated. The effect of the surface and 
bulk transverse fields on the phase diagrams of a semi infinite spin-1 ferromagnetic Ising model with a crystal 
field was investigated [118-119] within a finite cluster approximation with an expansion technique for cluster 
identities of spin-1 localized spin systems. 

Substantial research efforts have been devoted to the fabrication of quantum-dot structures, due to the in- 
teresting electronic and optical properties that can result from quantum conflnement. Several novel applications 
for such structures have been envisaged such as the single-electron transistor and intersubband quantum-dot 
infrared photodetectors. The well studied technique of creating quantum dots is the use of Si deposition on 
Ge [120,121]. A signiflcant number of studies have been performed with the aim of understanding the struc- 
tural and compositional evolution of Ge islands growing. The role played by different growth parameters such 
as temperature, pressure and others have been studied for different growth techniques. The evolution of self 
assembled Ge islands in low pressure chemical vapor deposition of Ge on Si, using high growth rates, has been 
investigated by atomic force microscopy and Rutherford backscattering spectrometry [122]. 

On the other hand, the understanding of the nature and consequences of adsorbate-adsorbate interactions 
at solid surfaces is of great interest in engineering properties of interfacial materials. Many recent studies show 
that indirect interactions, mediated by the structure, can be signiflcant enough to influence the formation of 
nano-structures at surfaces [123]. Indeed, much progress has been made recently in understanding adatom-pair 
interactions mediated by Shockley surface state electrons and the oscillatory interaction decays with adsorbate 
separation thickness [124-127]. Monte Carlo Simulations [128-130] and subsequently derived mean field theories 
[129,131], showed the existence of higher island densities than those expected by standard nucleation theory. It 
is shown that several experimental works, e.g. [127,132-134], supports these predictions, such as the deposition 
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of Au on mica substrates at a high temperature [135-142]. 

In most cases, the electrical conductivity and practical utility of discontinuous metal films, have been ex- 
tensively explored [143-145]. On the other hand, media that are magnetized perpendicular to the plane of 
the film, such as ultra thin cobalt films or multi layered structures [146,147], are more stable against ther- 
mal self-erasure [144,148] than conventional memory devices. In this context, magneto-optical memories seem 
particularly promising for ultrahigh-density recording on portable disks. The roughness and mobility of the 
magnetic domain walls [149,150] prevents closer packing of the magnetic bits, and therefore presents a challenge 
to reaching even higher bit densities. Increasing information storage to high densities may evolve through ex- 
tensions of current magnetic recording technologies. But such increases in storage density might be achieved 
by using other techniques such as holography, or micro machined nano-cantilever arrays [151]. The bit-writing 
with local probes may be thermally assisted by a current [152] or a laser beam that raises local temperature to 
the vicinity of the Curie temperature, resulting in the formation of a reversed domain with a rough wall. On 
the other hand, a strain implemented by a linear defect, in a magnetic thin film [153], can realize smooth and 
stable domain walls, that can be implemented without nano-scale patterning. Furthermore, composite materials 
such as multilayer coatings and isotropic nanocomposite coatings, having structures in the nanometer scale, can 
even show properties, which can not be obtained by a single coating material alone [154]. The most widely used 
coatings are TiN, TiC, TiCN, and combination thereof, as well as some coatings with lubrificating properties 
such as diamond-like carbon. Monte Carlo simulations are applied to study the equilibrium magnetic domain 
structure of growing ultra thin ferromagnetic films with a realistic atomic structure [155]. Near the percolation 
threshold a metastable magnetic domain structure is obtained with an average domain area ranging between 
the area of individual magnetic islands and the area of the large domains observed for thicker ferromagnetic 
films. This micro-domain structure is controlled and stabilized by the non-uniform atomic nano-structure of 
the ultra-thin film, causing a random interaction between magnetic islands with varying sizes and shapes. The 
investigation of the magnetic domain formation, in nano-structured ultra-thin films, is an active field of current 
research. The influence of the atomic morphology on the magnetic properties is known to be especially strong 
during the initial states of the thin film growth. A small variation of the preparation conditions may consid- 
erably change the obtained magnetic structure. This has been shown, by resolving imaging techniques, which 
allows the investigation of atomic structure as well as the magnetic domain structure of several nanostructured 
systems [156,160-161]. The influence of the presence of a surface up on the bulk critical behavior has been the 
focus of scientific interest during the last decade [162-164]. Since the free surface is the natural break of the 
translational symmetry in any real system, critical phenomena at perfect surfaces have attracted much interest. 
Furthermore, the phase diagram for the semi-infinite Ising model is well established [162,164]. Moreover, surface 
critical exponents have been estimated, both theoretically. The general agreement is quite satisfactory [162-164]. 
The critical behavior of magnetic materials with variable strength of coupling on the surface is characterized by 
surface critical exponents independent of the bulk ones. For their theoretical determination e-expansion [164], 
series expansion [165], and Monte-Carlo simulations [166,167] were applied. Moreover, the magnetic properties 
depend on several geometric features, which depend on the conditions during the growth, thus, the connection 
between experimental and theoretical results requires a detailed study of the dependence of physical quantities 
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on the geometry of the systems. Some specific properties have aheady been considered in theoretical works, 
such as the size and the lattice structure of small clusters [168,169], the thickness of the film [170] and the 
roughness of one dimensional structures [171], in particular, the role of surface imperfections, unavoidable in 
real materials, on the critical properties has been studied [172-175]. According to the Harris criterion [176], the 
relevance of the perturbationis connected to the sign of the specific heat exponent a in the pure system. Small 
amounts of disorder do not change the nature of the phase transition,if the corresponding pure system's specific 
heat exponent is negative. Currently, there appears to be widespread consensus that critical exponents of bond 
disordered systems in the weak disorder limit are the same as those of the pure one [177-179]. Monte-Carlo 
simulations of bond disordered systems show logarithmic corrections in the behavior of the magnetizations 
which persist even in the strongly disordered regime [180,181]. However, in the 2D spin diluted Ising model a 
continuous variation of critical exponents with the spin density has been found [182]. A Monte-Carlo study [183] 
shows that the specific heat does not diverge at the critical point, at least when the concentration of dilution is 
sufficiently large; the critical behavior of the magnetic susceptibility and the correlation length are consistent 
with the pure power-law behavior with the concentration dependent critical exponents. Moreover , using Monte 
Carlo simulations, the site-diluted Ising model in two dimensions has been studied [184] where it is found that 
such model belongs to the same universality class as the pure one, in spite of the strong logarithmic efii'ects. 
Besides, the surface critical behavior of the three-dimensional disordered semi-infinite Ising ferromagnetic has 
been studied, using Monte Carlo simulations [185]; it is found that for less than Jgc, the corresponding surface 
critical exponent is equal to the one of perfect surface [186-192]. 

Using the perturbative theory, Harris [193] have studied the layering transitions at T = in the presence of 
a transverse field. The effect of finite size on such transitions has been studied, in thin film confined between 
parallel plates or walls, by Nakanishi and Fisher [194] using mean field theory and by Bruno etal. [195] taking into 
account the capillary condensation effect. The critical wetting transitions have been considered in systems with 
corrugated periodic walls [195] by using an effective Hamiltonian study. By applying Monte Carlo simulations 
on thin Ising films with competing walls, Binder et al. [196], found that occurring phase transitions belong to the 
universality class of the two-dimensional Ising model and found that the transition is shifted to a temperature 
just below the wetting transition of a semi-infinite fiuid [197,198]. Hanke etal. [199] show that symmetry breaking 
fields give rise to nontrivial and long-ranged order parameter profiles for critical systems such as fluids, alloys, 
or magnets conflned to wedges. The existence of the layering transitions for a film, with corrugated surface 
has been studied using mean field theory and Monte Carlo simulations [200]. It is found that the wetting 
temperature Tyj depends weakly on the surface corrugation degree for fixed values of the surface magnetic field. 
In addition, by applying a suitable effective interface model at liquid-vapor coexistence, Rejmer etal. [201] found 
a filling transition at which the height of the meniscus becomes macroscopically large while the planar walls of 
the wedge far away from its remain non-wet up to the wetting transition occurring at T^. They also showed 
that the discontinuous filling transition is accompanied by a pre-filling line extending into the vapor phase of 
the bulk phase diagram and describing a transition from a small to a large, but finite, meniscus height. Much 
attention has been paid to the properties of layered structures consisting of alternating magnetic materials. 
The most commonly studied magnetic multilayers are those of ferromagnetic transition metal such as Co or 
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Ni. Many experiments have shown that the magnetization enhancement exists in multilayered films consisting 
of magnetic layers. It was found that ferromagnetic coupling can exist between magnetic layers. Prom the 
theoretical point of view, great interest has been paid to spin wave excitations as well as critical phenomena 
[202-215]. The study of thin films is partly motivated by the development of new growth and characterization 
techniques, but perhaps more so by the discovery of many exciting new properties, some quite unanticipated. 
These include, more recently, the discovery of enormous values of magnetoresistance in magnetic multilayers far 
exceeding those found in single layer films and the discovery of oscillatory interlayer coupling in transition metal 
multilayers. These experimental studies have motivated much theoretical work. However these developments 
are applied to a large extent powered by "materials engineering" and the ability to control and understand 
the growth of thin layers only a few atoms thick. However, Many experimental studies have shown that the 
magnetization enhancement exists in multilayered films consisting of magnetic layers. The development of the 
molecular beam epitaxy technique and its application to the growth of thin metallic films has simulated renewed 
interest in thin film magnetism. The fabrication of thin metallic films has led to a surge of experimental activity 
[216-219]. More recently, the thickness dependence of the critical temperature of thin Ising magnetic films 
having surface exchange enhancement has been studied by Hong [71] and by Aguilera-Granja and Moran Lopez 
[220] within mean field approximation, by Hai and Li [221] and Wiatrowski et al. [222] within an effective field 
treatment that accounts for the self-spin correlations. 

On the other hand, ferroelectric films can be described by an Ising model in a transverse field [223-225], 
and when the film becomes very thick, its properties are those of the semi-infinite Ising system [108-110,226]. 
Phase transitions in Ising spin systems, driven entirely by quantum fluctuations, have been getting a lot of 
attention [227]. The simplest of such systems is the Ising model in a transverse field which can be exactly 
solved in one dimension. This model has been introduced to explain the phase transition of hydrogen-bonded 
ferroelectrics such as KH2PO4 [228]. Since then, this model has been applied to several physical systems like 
DyV02 and studied by a variety of sophisticated techniques [92,96,229,230] .The technique of differential operator 
introduced by Kaneyoshi [231] is as simple as the mean field method and uses a generalized but approximate 
Callen relation derived by Sa Barreto and Pittipaldi [115]. The system has a finite transition temperature, which 
can be decreased by increasing the transverse field to a critical value Oc- The effect of a transverse field on 
the critical behavior and the magnetization curves was studied [92,96,229-233]. Using the perturbative theory, 
Harris et al. [193] have studied the layering transitions at T = in the presence of a transverse field. The 
effect of a uniform transverse magnetic field on the layering and wetting transitions of a spin— 1/2 Ising model 
in a longitudinal magnetic field and showed the existence of layering and wetting transitions above a critical 
transverse field Q.^, which is a function of the temperature and surface magnetic field [79] . Recently, Karevski 
et al. [82] have studied the random transverse Ising spin chain according to a distribution of the transverse field 
governed by a law of type: z~"; z being the distance from the surface and a a constant). 

On the other hand. The simplest of all random quantum systems is the random transverse Ising model 
[234,235]. However, the random systems have been known to be dominated by rare regions. This effect is 
particularly pronounced for random quantum systems at low or zero temperature far from critical points. 
Indeed, Griffiths [236] showed that the free energy is a non analytic function, because of rare regions. Having 
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found all the derivatives being finite, Harris [237] concluded that this effect was very weak for classical systems. 

The spin-1 Ising model with nearest-neighbor interactions, both bilinear and bi-quadratic, and with a crystal- 
field interaction was introduced by Blume, Emery, and Griffiths (BEG) [104] to describe phase separation and 
superfluid ordering in ^He—^He mixtures. With vanishing bi-quadratic interactions the model is known as the 
Blume-Capel model [103]. Since its introduction, the spin-1 BEG model has been extended to solid-liquid-gas 
systems, and multi-component fluid liquid crystal mixtures [105], magnetic materials [238-240], critical behavior 
and multi-critical phase diagrams [241-244]. Furthermore it is the simplest model that can be used for modeling 
the behavior of the liquid, solid and vapor phases of a real materials. In this context it was used by Jayanthi 
[245] to study the surface melting of the solid phase near the triple point and by Gelb [246] to characterize the 
layering transitions and surface melting both near and away from the bulk coexistence lines under the effect 
of temperature. The effects of quantum transverse field on wetting and layering transitions of a spin-1 BEG 
model are also studied [79,247,248]. 

2 Effects of surface potential and surface coupling on wetting and 
layering transitions 

2.1 The effect of magnetic surface field 
2.1.1 Model 

We consider a spin-1/2 Ising model, with ferromagnetic coupling on a three dimensional cubical lattice consisting 
of N parallel square lattice layers at spacing a to make a "film" of finite thickness , limited by a single substrate. 
Each layer has A sites, labeled i,j , and at each site is an Ising spin variable (Ji, which is in the simplest 5 = 1/2 
case takes the values ±1. In the simplest model, only nearest-neighbor interactions are considered and the 
Hamiltonian is then 

N N-1 N 

w = - E E -^^.i^^ - E •^^.^'+1 E ^''^f^' - E E (1) 

p=l p=l i p=l i 

where is the sum over layers, "^^j means a sum over distinct nearest-neighbors pairs on the layer p 
and J2i means a sum over nearest-neighbors sites between adjacent layers. While indicates a sum over all 
sites of the lattice. N is the number of layers. Hp, the magnetic field applied to the layer p. Jij = Jp for 
located in the layer p ( the intra-layer coupling), while Jp,p+i is the analogous interlayer coupling. In the 
uniform magnetic field case Hp is given by: 

Hp = HJp^i + H (2) 
with Jp^i = 1 for p = 1 and zero elsewhere. While in the variable magnetic field case Hp is given by 

Hp = ^+H (3) 
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Hs is the surface magnetic field, which depends on the nature of the substrate and H is the external magnetic 
field. For technical reason it is practical to consider the reduced Hamiltonian as follows: 

-m = -j2K,Y. <^>'o - E E -flKj:^! (4) 

p=l p=l i p=l i 

(3 = j^^, Kp = (3Jp and hp = (3Hp; Kb is the Boltzmann constant and T is the absolute temperature . 
2.1.2 Real space renormalization group technique 

The Hamiltonian in the Eq.(4) is not invariant under the Migdal-Kadanoff renormalization group transfor- 
mation. Then, it is more appropriate to redefine (3'H in terms of a sum of general 2A''-sites interactions: 
f^,j{alaf...af;a}al..af) so that 

m = -Y: ; (5) 

The sum is taken over nearest-neighbor sites on a square lattice for the three-dimensional cubical lattice, and 
over nearest-neighbor sites on a chain for the two dimensional square lattices. h^^j{ajaf...aj^;ajaj...aj'), 
where n is the renormalization-group transformation index given initially by : 

N N-1 N 

= - E - \ E <i+i(-r-r^ + <-r^) - - E <H-r + o (6) 

Where ^ = 4 in the three dimensional cubical lattice case, and ^; = 2 in the two-dimensional square lattice 
.Kf^ = -f^p"p_|_i = PJ and hp'^ = hp are the parameters of the system before renormalization. Eq.(6) has 
the simplicity of formally having reduced the N-layer system to an effective one layer system with a "2N-site" 
nearest-neighbor coupling. Hence, the partition function is given by: 

Z = Tr,^ H exp{-hfl,.) (7) 

We now invoke a Migdal-Kadanoff bond-moving transformation in order to reduce systematically the degrees 
of freedom in Eq. (7). We may use the displacement of bonds on square lattice, with a scale factor b=2. Then 
we sum over the variables associated to the crossed sites[69,211] (see Fig. 1). The recursion relation is easily 
found to be : 

m >{(T^(T^...a^ :a^cT^...a^ \ (8) 

with e = 1 in the two-dimensional square lattice case and e = 2 in the three dimensional cubical lattice case; 
is the normalization constant, where we have used the reflection symmetry, 

'{aiai...ai ■,ajaj...aj ) = U"- '{a^aj-.a^ \<Jiai...<Ji ) (9) 

Note that Eq.(8) provides one with a set of 2^^ recursion relations since each variable (Ji can takes two 
values ±1. The reflection symmetry (9) reduces this number to 2'^^~^ + 2^~^ , one of which is the normaliza- 
tion condition. In order to extract useful information from iterations of Eq.(7), we write the general 2N-site 
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Figure 1: The Migdal-Kadanoff Model for 6=2. 



Hamiltonian in the form 

2JV /(r) 

- ...af ; a]oj...af) = ^ i^^ (10) 

r—l m=l 

where f(r) is the number of different r-site Ising interactions, and 

2JV 

^/(r) =22^-1 + 2^-1-1 (11) 

r=l 

(n) 

Km'r is the interaction of type m between r spins, is the function of type m of r coupled spins, which is 
symmetric by the permutation of i and j. As an example, we give the simplest case N = 2, where 



- }^j{ala^;a]aj) = Ki^h^l + a]) + {a^ + a]) + Ki^^ {aj aj) 
+/^g(a?aj) + K^ialaf + a]a]) + K^:i{a\a] + a^af) 
+i^g(aiaf a| + + Ktl{a\a]a] + a]a\a1) + K^rlo\a]ala] (12) 

with k'^^I and k'^}^ are respectively the rcnormalizcd reduced magnetic field; k'^i^ are respectively 

in] 

the intra-laycr rcnormalizcd reduced coupling in the first and the second layer; K\ 2 is the intcrlayer rcnor- 
malizcd reduced coupling ■.k'^^^ t ^'Tz ^ ^23 '^^^'^ ^i^I (hi^itially equal to zero) arc the rcnormalizcd reduced 
coupling between two, three and four spins, which permit to have an homogeneous system of equations. Our 
rcnormalization group calculation show in all situations that the reduced intcrlayer coupling iterates to infin- 
ity for the two-dimensional layers, while these parameters iterate to a fixed point (at T = 0, this parameters 
iterate to infinity) in the one-dimcnsional layer case. The behavior of the other parameters in Eq. (10) is 
much more interesting in that it reflects the transition foimd in the three-dimensional ferromagnet films [69,70]. 
The ferromagnetic layer state is characterized by the value of the rcnormalizcd reduced magnetic field, applied 
on this layer, in the trivial fixed point. If Zip"' iterates to +00 for sufficiently large iteration n, the state of 

1. . (n) 

the layer p is ferromagnetic "up" and if hp iterates to —00, the state of the layer p is ferromagnetic "down". 
We remark that the result is independent of the choice of the scale factor of the real-space rcnormalization group. 



2.1.3 Transfer matrix method (TM) 

The model described by the Hamiltonian (1) is also studied using transfer matrix method . This method 
consists in studying a d-dimensional strips which are infinite in one direction and finite in the remaining ones 
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Figure 2: The strip used within the finite size method, the transfer is made from the plane I to the plane l + l. 
Cyclic boundary conditions are made in the direction k and free boundary conditions in the direction p (p is 
the layer index). 



[70] . This system with a periodic boundary condition on one finite direction (parallel to the substrate) and the 
free boundary conditions on the perpendicular direction to the substrate, is similar to a (rf — 1) dimensional 
layers (see Fig. 2). The partition function Z of this strip is written under the matrix product form (2^^, 2^^) 
called the transfer matrix, M is the width of the strips and N is the number of layers, namely: 

oc 

Z=l[V{L,L+l) (13) 
This makes the transfer, from the plane L to the plane L+l, in the infinite direction. 

The elements of this matrix are given by 

V{i,j) = exp(^a{p,k,i){a{p + l,k,i) + a{p, k + + a{p, k, j)) +^hpa{p,k,i)]/t) (14) 

p,k p 

where hp = Hp/ J; t = T/J and a{p, k, i) is the spin variable which takes two values ±1; p = 1, N; i and j 
takes values ; p is the layer index which takes N values and k takes M values. The periodic conditions 
on the direction k indicate that 0-(A'' + l,k,i) = 0. The free energy per site is given by / = —TLog{\Max) 
where Amox is the largest eigenvalue of the transfer matrix. In order to determine the state of each layer, it is 
necessary to define the magnetization per site on the layer p as follows, 

where hp is the reduced magnetic field applied on a layer p. The magnetization nip of a layer p can be positive 
or negative. However, if rup is positive, the state of the layer p is ferromagnetic "up", while if vrip is negative, 
the state of such layer is ferromagnetic " down" . The transition of this layer is given by the change of the 
magnetization sign. 
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2.1.4 Mean Field Theory 

The mean field state equations of the Hamiltonian (1) may be obtained from the variational functional, 
F[mp] = T[-^Log{l - nip) + -Ip^Log{l + rrip) - Log{2)] 

-Hpmp - ^mp{mp-i + Amp + nip+x} (16) 

with the free boundaries conditions toq = rriM+i = 0. At the global minimum of the functional free en- 
ergy(Eq.l6), the magnetization per site in the layer p is < >= m™*", and the free energy of the system is 
F = F[rn™'"]. The stationary condition dF/dnip = leads to the mean field state equations, 

mp = tanh[^{mp-i + ArUp + mp+i + hp)];p = 1, N (17) 

These equations are solved numerically using iterative method. Diflferent initial guesses can lead to different 
solutions of the Eq. 17; the one which makes the variational functional free energy smallest is selected as the 
global minimum. At the transition of the layer p, the magnetization rUp has a discontinuity for a temperature 
smaller than the critical-end-point temperature. 

2.1.5 Uniform surface magnetic field 

At first it is useful and instructive to find the ground state of the Hamiltonian (1) describing the phase diagrams 
in the {hg — h) plane at T = QK, this is determined analytically. We denote here after by l^'O^ the configuration 
having p layers with positive magnetization and q layers with negative ones. The ground state phase diagram 
is presented in Fig. 3. It is found that there exist two different situations depending on the value of the surface 
magnetic field. Indeed, for hg < N/{N — 1) we have one transition (0^ 1^), a,t H = —hg/N in which the 
N layers transit simultaneously and for hg > N/{N — 1), we have two transitions; one a,t H = —1/{N — 1), 
(1^ <-> 10^~^) where the surface magnetization does not changes the sign, but the remaining {N — 1) layers 
transit simultaneously , and the second (10^~^ 0^),at H = a + bhg, where only the top layer (surface) 
transits (i.e the surface magnetization changes the sign). 

At finite temperature, numerical calculations are done in the case of a thin films N = 3. However, within 
RSRG and TM methods we obtain one transition (0^ ^ 1^) for hg < 1.5 and T <Tl (Fig. 4a), in agreement 
with the ground state results, and three transitions (0^ 10^ 1^0 ^ 1^) for Tl < T < Tc. For hg > 1.5, 
we obtain two transitions (0^ 10^ 1^) for T < Tl as predicted by the ground state calculations and 
three transitions (0^ 10^ ^ 1^0 1^) (Fig. 4b). Tc is the transition temperature of a three layer film 
(A'' = 3). The layering temperature Tl above which a sequence of layer transitions occurs depends strongly 
on the strength of the potential of substrate, i.e the value of the surface magnetic field hg. Within MFT only 
the ground state phases appear at finite temperature, since MFT neglects correlation between spins, hence the 
effects of thermal fluctuations are negligible, and then does not raise the degeneration of the states at T > 0. 
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Figure 4: Phase diagram in (T — ft,) plane for iV = 3 obtained within RSRG in the uniform magnetic case ; (a) 
hs = 1.2 < 1.5 ; (b) hs = 1.8 > 1.5 [69,70] 
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2.1.6 Variable surface magnetic field 

In this case the problem is quite different from the uniform surface magnetic field since for sufficiently large 
value of hs the layering transitions occur at T = OK (Fig. 5). However for hsg < hg < hsq+i with q= 1, ...TV ; 
we have q transitions in the hs — h plane namely (0-^ 10-^"^ 1^0^~^ ...l'~^0^~'^ 1^), where hs-^ = 0, 
= 1.63, = 11.36 and hg^ = oo. 

Using RSRG and TM methods, three different cases are distinguished, at finite temperature: 

i) For hs < hs2 only one transition (0^ 1^) persists for T <Tl and three transitions (0"^ 10^ ^ 1^0 1^) 
for Tl<T <Tc (Fig. 6a). 

ii) For hs2 < hs < hs^ we have two transitions (0^ ^ 10^ ^ 1^) for T < Tl, while for T > Tl, (Fig. 6b) we 
obtain three transitions (0"^ ^ 10^ ^ 1^0 <-^- 1^). While MFT gives rise to only two transitions (0^ 10^ 1^) 
for any value of hs smaller than hg-^ even for T > Tl- 

iii) For a strong potential of the substrate {hs > hg^), the layering temperature Tl vanishes, and then, the 
three methods give rise to the appearance of three transitions, (0^ ^ 10^ <-» 1^0 <-» 1^), for any temperature 
0<T<Tc (Fig. 6c) 

2.2 Surface coupling effects 

We consider the Hamiltonian given in Eq.l, in which, in this case, the surface and bulk couplings are different; 
Ji,j = Js and Jij = Jh if {hj} are in the surface and the bulk respectively (Fig. 7). 

In the following, the notation l''0^~'' means that a configuration of the system have the first k layers with 
positive magnetizations and the remaining N—k layers with negative ones. However, there exist A'' + 1 possible 
configurations for a film formed with A'' square layers. Numerical calculations are limited to a = 2.0, since the 
results for other values of a are qualitatively the same. 

At T = 0, and for {N > 3), there is only three allowed transitions. ^ 1^, ^ lO^'^ and lO^'^ ^ 1^. 
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Figure 6: Phase diagram in (T — h) plane for A'' = 3 in the variable case ; (a) hs = 1.2 ; (b) /is = 4 ; (c) hs = 12. 
Sohd Hnes correspond to the RSRG method and dashed Unes correspond to MFT [69] 
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Figure 7: A sketch of the geometry of the system formed with layers and subject to the surface magnetic field 
Hs and the global magnetic field H . The distribution of the coupling constants : Js denoting the interaction 
coupling constant between the spins of the surface and J5 corresponding to the interaction coupling constant 
between the spins of the bulk as well as between the spins of the bulk and those of the surface ; are also 
presented. 
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Figure 8: Phase diagram, in the {H/Ji,,T/Jh) plane, showing the layering transitions, and the definition of the 
layering temperature T^/Jb, which is different from the surface transition. This figure is plotted for a system 
size with = 10 layers, a magnetic surface field Hs/Jb = 1-0 and a surface coupling constant Js/Jb = 1-5 [250] 



At T 7^ 0, Fig. 8 shows that the layering temperature, Tl, above which a succession of layering transitions 
occurs is greater than the "surface transition" temperature. Hence, the film thickness must be at least equal 
to three layers and so the results established throughout this work correspond to films with a thickness N > 3 
layers. 

It is understood that the wetting temperature T^/Jb is close to the limit of TL/Jb for a sufficiently thick 
film {N — > oo). The layering temperature behavior as a function of the system size is illustrated in Figs. 9a 
and 96 for two surface field values and several values of the surface coupling. However, for smaller values of the 
surface field, the layering and wetting temperatures are independent of the surface coupling for large system 
sizes as it is shown in Fig. 9a. While for higher values of the surface field (Fig. 9b), the wetting and layering 
temperatures depend strongly on the surface coupling values. Indeed, there exist three different classes of the 
layering temperature behaviors: 

(i)Tj, increases with the film thickness and stabilizes at a certain fixed value. Such behavior occurs at small 
values of Jg and any value of hg. 

{ii)TL increases until a certain film thickness above which it decreases exhibiting a maximum. This situation 
take place for medium values of both Jg and hs- 

{iii)TL decreases continuously and stabilizes at a value close to T^. This behavior occurs at higher values of Jg. 

Furthermore, the wetting temperature T^, which is the limit value of Tl for L sufficiently large, is indepen- 
dent of Js for small values of hs- While for higher values of hs, it depends strongly on the Jg and decreases when 
hs increases. Indeed, this is because the surface region becomes magnetically harder than the deeper layers, 
and as the film is thicker, the strength of the internal magnetic field decreases, which leads to a decrease of T^. 
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Figure 9: Layering temperature, T^/Jf), behavior as a function of the number of layers, a) Hgj J\, = 0.1 and b) 
Hg/Jf) = 1.0. The number accompanying each curve denotes the surface coupHng Jg/J}, value [250] 



We note that, the three distinct behaviors mentioned above are separated by two critical surface coupling J ad 
and Jsc2- Indeed, J ad separates the behaviors (i) and (m), while Jac2 separates the behaviors (ii) and {iii). 
The dependence of Jgd and Jac2 as a function of hg is plotted in Fig. 10. 

Jsd and Jsc2 correspond to two critical surface fields had and hsc2 separating the three different regimes 
of the layering temperature T^. However, for sufficiently large values of J, ( Js/Jb > 1.30), decreases when 
increasing the film thickness N for any values of ha. While, for Jg < 1.3, the temperature Tl exhibits the three 
behaviors; (i), (ii) and (iii) depending on the value of kg. In particular for ha > hac2, is always decreasing 
function of the film thickness. The layering temperature behavior is plotted, in Fig. 11, for a fixed system size 
N = 5 layers and several value of ha . It is found that increases and decreases for small and higher values of 
hs respectively. On the other hand, the wetting temperature, Fig. 11, is not affected by the surface coupling 
variations for small surface field values, while for higher hs, decreases abruptly and stabilizes at a certain 
limit, when increasing the surface coupling. 

2.3 The surface crystal field effects 

The system we are studying here is formed with N coupled ferromagnetic square layers in the presence of a 
crystal field. The Hamiltonian governing this system is given by 

w = - E -^^i^-^o + E ^*(^^)' (18) 

<i,j> i 

where, Si{l = = —1, 0, +1 are the spin variables. The interactions between different spins are assumed to be 
constant so that Jij = J. All sites i of a given layer k, {k = 1, 2, N) are subject to a crystal field amplitude 
Aj = Afc so that Ai > A2 > ... > Ajv- The crystal field of a layer k is given by: 

Afc = A,/A:" (19) 

where A^ = Ai is the crystal field acting on the surface (first layer k = 1) and a a positive constant. For a system 
with a = 0, a constant crystal field is applied on each layer of the film. Whereas, for a — > 00, the crystal field 
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Figure 11: Layering temperature and wetting temperature beliaviors, for the (USF) case, as a funetion of the 
surface coupling Jg/Jt and selected values of the surface magnetic field iJ^/Jf,: 0.1, = 0.2, 0.3, = 0.4, 0.5 and 
0.6. Ti/Jb is plotted for a fixed system size with N = 5 layers [250] 
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is applied only on the first layer (fc = 1). Experimental studies of a sequence of reentrant layering transitions of 
argon and krypton adsorbed on graphite substrate are given in Ref . [54] , in which the algebraically decaying of 
the applied field according to a law of type (A:~")have been used. However, the quantities computed, using the 
Monte Carlo simulations for each layer k containing spins in a;— direction and A^^ spins in the ?/— direction, 
are: 



-The layer average magnetizations 



-The layer magnetic susceptibilities 



mk = iJ2s^)/iN,Ny) (20) 

iek 



Xm,k = PN^Ny < {rrik- < nik > f > (21) 
-The layer quadrupolar magnetic susceptibilities expressed as 

Xg.fc = PN^Ny < {qk- < Qk >)' > (22) 

with 

qk = {^Si^)/{N,Ny). (23) 

iek 

-The layer critical exponents 'ym,k related to the corresponding layer magnetic susceptibilities Xm,k, for a fixed 
layer k, giving by 

dLog\T,,k/J-T/J\ ^^^> 
where T/J stands for the reduced absolute temperature and T^^k/J is the reduced critical temperature of the 
layer k. In the above equations /3 = l/{kBT) with ks being the Boltzmann constant. 

The notation D^O^'^ means that the first k layers from the surface are disordered while the remaining N ~ k 
layers are ordered. The ground state phase diagram (Fig. 12) of this model [102], shows that, there exist only 
fco possible layer transition which depends on a and N. 

For T > Monte Carlo simulations shows, for N = 10, that the last layers fc = 8, 9 and 10 transit 
simultaneously at As/ J « 23.84. 

Phase diagrams exhibit critical, tricritical and reentrant behaviors (Fig. 13). The latter is due, for the first 
ko layers, to the competition between thermal fluctuations, the crystal field and an effective magnetic field 
created by the deeper ordered layers. Indeed, when these thermal fluctuations become sufficiently important, 
the magnetization of some spins, of deeper layers, becomes nonzero (+1 or -1). This leads to the appearance of 
an effective magnetic flcld. 

This magnetic field is responsible of the ordered phase seen for the layer k. This argument can also explain 
the absence of the reentrant phenomena for the last N — ko layers, once the value of the crystal field is not 
sufficient to overcome the effective magnetic field created by the remaining deeper layers. It is worth to note 
that the reentrant phenomena is always present for the layers fc, (fc < fco), and the corresponding tri-critical 
points Ci are located at a constant temperature. At the first order transition the layer quadrupolar magnetic 
susceptibilities Xq.k , present a strong peak at the surface crystal field transition (Fig. 14). 

While the magnetic susceptibilities Xm,kj present a strong peak at the second order transition (Fig. 15). 
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Figure 13: The critical temperature behavior as a function of the surface crystal field Ag/J for a = 1.0 and a 
film thickness A'' = 10 layers. For each layer k {k = 1,2, ...,N), the first-order transition line (vertical dashed 
line) is connected to the second-order transition line (up-triangular points) by a tri-critical point Cfe (open 
circle) [251] 
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Figure 14: The behavior of the layer quadrupolar magnetic susceptibilities Xq,k, showing a first-order tran- 
sition, as a function of the surface crystal field Ag/J, for a very low temperature: T/J = 0.2. The number 
accompanying each curve denotes the layer order from the surface to deeper layers [251] 
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Figure 16: Thermal behavior of the layer magnetization my (line with solid circles) and the reduced layer 
magnetic susceptibility Xm,7 (line with open circles) for a deeper layer k = 7. for As/ J = 22.0 [251]. 



In the case of reentrant behavior, Xm.k present two successive peaks at the transition points (Fig. 16). The 
critical exponents jm.k of magnetic susceptibility of each layer, for N = 10. 

It is found that jrn.k decreases, for a fixed order fc, with the system size x Ny and stabilizes at certain 
value; and due to the free boundary conditions the critical exponents of the layers k ~ 1 and k = N arc found 
to be greater than those of the internal layers 2 < fc < iV — 1. We note that the critical exponent of a 2D 
layer in a film is different than the one of a free 215-layer. Indeed the critical exponent of a layer depends on its 
position (Fig. 17). 

2.4 The effect of the inhomogeneity of substrate 
2.4.1 Model 

The effect of the inhomogeneity of the substrate on the wetting transitions is studied, using Monte Carlo 
simulations, in the case of a spin-1/2 Ising ferromagnetic film formed by N coupled square layers, in which the 
surface magnetic field acts only on alternate clusters of spins of the surface (fc = 1) indicated by the symbols 
(+); while Hg is absent on clusters of symbols (o) (Fig. 18). The Hamiltonian governing this system is given by 

H = -JY, S,Sj - + Hs,)Si (25) 

<i.J> i 

where, Si{l — = — 1,+1 are the spin variables. The interaction between different spins is assumed to be 
constant. H is the external magnetic field. 

The surface magnetic field if^. applied on each site i of the surface A; = 1, is distributed alternatively, so 
that: 

{+Hs for all sites i e clusters with symbols (+) 
(26) 
for all sites i e clusters with symbols (o). 
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Figure 18: A sketch of"tTie°syst^m^edmetry for A'' = 4 layers. A surface magnetic field Hs is acting on 
alternate island spins of the surface k = 1. Hg is present on clusters with symbols: '+', and absent in clusters 
with symbols 'o'. 
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Figure 19: The (T, H) phase diagram for the layer fc = 4, for Hg: 1.5; 2.0 and 2.5 [252]. 



2.4.2 Phase diagram and separation of holes transitions 

We mean by hole a cluster formed by positive spin surrounded by negative one. However, depending on the 
values of H, the (T, H) phase diagram established in Fig. 19 shows the existence of three different phases, for 
a fixed surface magnetic field Hg. The totally wet (TW), the nonwct (NW) and the partially wet phases. 

In the latter case, we can distinguish, in each layer, three different holes configurations namely; the (PWTD) 
configurations in which the layer is partially wet with a total disconnection of the holes; the (PWPD) configu- 
rations in which the layer is partially wet with a partial disconnection of the holes; the (PWTC) configurations 
corresponding to a partially wet layer with a totally connected holes. Fig. 19 shows that at sufficiently low 
temperature, the Partial wet disappears, while at higher temperature values, provided that T is kept less than 
the critical value Tc = 3.87, the partially wet phase is found. 

Furthermore the surface magnetic field Hg favors the appearance of the (PWPD) phase. In Figs. 20, we 
give, for the fourth layer, an example of three different partial wet phases; (PWTC), (PWPD) and (PWTD) in 
which the distributions of holes are different. 

2.4.3 The size distribution of holes in a layer 

In the (TW), (PWTC) and (NW) phases, all spins are negative, only one big hole exist and all spins are positive 
respectively. However the distribution of size of holes is important in both (PWTD) and (PWPD) phases in 
which there is a great number of disconnected holes with different sizes. In Fig. 21 we give the distribution 
of sizes of holes, in each layer {k = 1,2,3,4) of the film, in the (PWTD) phase, for a temperature T = 3.5, 
Hg = 2.0 and H = —0.142. It is clear that, there exist three different peaks at three different sizes 2, 40 and 
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Figure 21: Frequency islands size distributions for the layers k = 1,A; = 2,A; = 3 and fc = 4 for the configuration 
(PWTD) for T = 3.5, = 2.0 and H = -0.142 [252]. 

60. The difference between the size frequencies decreases when increasing the value of k. Indeed, far from the 
surface, the effect of the surface field becomes less relevant and then the holes are equally distributed. Such 
result is illustrated in Fig. 21 in the A; = 4 case. 

3 Geometry effect on wetting and layering transitions 

3.1 Corrugated surface effects 
3.1.1 Model and method 

We consider a fihn with finite thickness of a three dimensional spin-1/2 Ising model, with M square layers, 
limited by two surfaces (Fig. 22). The top one is corrugated with n steps, each step contains L spins. The top 
corrugated surface St is under a uniform magnetic field Hst and the bottom surface Sb is a perfect plane with 
an applied uniform magnetic surface field Hst- 

The Hamiltonian of the system, can be written as 

H = — JijSiSj — ''^^HiSi (27) 

<i,j> i 

where. Si = ±1 are spin— 1/2 Ising random variables, and Jij = J arc the exchange interactions assumed to be 
constant. The system is assumed to be infinite in the direction y, so the variable y will be cancelled in all the 
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Figure 22: Geometry of a corrugated surface system with M layers in the 2;— direction and n steps, each step 
contains L spins in the a;— direction. The system is infinite in the y— direction. A surface field Hst is applied 
on the top surface St, whereas the bottom surface Sb is under a surface field Hst = —Hst- An external field H 
is applied to the system. 
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Hi = H{x, z) = < 



following. Hi is the total longitudinal field applied on a spin located on a site 'i' of coordinates {x,z). defined 
by: 

H + Hst for ieSt 

H for ieBulk (28) 

H + Hsb for ieSb 

H is the external magnetic field. In the following we will limit our interest to a system with Hsb = —Hst- 
The mean field equations describing the magnetization of the system is given by 

m{x, z) = tanh(/3(2m(a;, z) + m{x, 2; + 1) + m{x, 2; — 1) + m{x + 1, z) + m{x —\,z) + H(x, z))) (29) 

The total free energy of the system can be derived as : 

F[m{x, z)] = J2, jT['-^^Log{l - m(x, z)) + ^-±2Ip^Log{\ + m(x, z))\ 

~jm{x, z)[2m{x, z) + m{x, z + 1) + m{x, z ~ I) + m{x + 1,2;)+ m{x — 1, z)] (30) 

—m{x, z)H{x, 2)} 

The free boundary conditions are considered to solve the mean field equations. 
3.1.2 The corrugated surface critical behavior 

In this case the critical exponents are computed, for iJg = and H = 0, using mean field theory. At the vicinity 
of the critical temperature Tc, the spontaneous magnetization m{z,t), of each layer z, vanishes with a power 
law given by 

m{z,t) =t^-ft'^^'*^ (31) 
Then the effective exponent f3eff{z,t) can be defined as 

I3eff{z, t) = dln{m{z, t))/dln{t) (32) 

where t = \T — T^jT^ is the reduced temperature and m{z,t) is the average value of ■m{x,z,t) on x. For 
sufficiently small values of t, f3effiz,t) approaches the asymptotic critical exponent (i{z). The surface critical 
occurs when 2; = 1 denoted hy pi = P{z — 1). The effective exponent can be approximated at discrete reduced 
temperatures ti and as 

f3effiz, t) = ln[m{z, ti)/m{z, ti+i)]/ln{ti/ti+i). (33) 

However, in the vicinity of the film critical temperature T^, the profile of the effective critical exponent f3eff{z, n) 
is presented in Fig. 23. 

It is clear that at the surface, the cfFoctivc critical exponent remains independent of the number of steps. 
While for any layer z with ni < 2 < n2, it depends strongly on the number of steps. For a fixed number of 
steps, the thermal variation of the effective exponent is presented in Fig. 24. 

It is found that far from the transition temperature, (3eff{z, n) decreases with increasing z, shows a crossover 
&t z — n, continues to decrease until it becomes constant in the depth of the system and increases again in the 
vicinity of the bottom surface. If we increase the temperature to approach the critical temperature, the crossover 
becomes less pronounced. The dependence of the effective critical exponent as a function of temperature shows 
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Figure 24: The dependence of the effective exponent as a function of reduced temperature t, for L = 3f , d = 60, 
n = 4, and Ji = Js = Jl- The number accompanying each curve denotes the value of the layer position z [198] 
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Figure 25: The dependence of the effective exponent as a function of the number of steps n, for L = 31, d = 60, 
T = 0.90Tc and J\ = Jg = Jl- The number accompanying each curve denotes the layer position z [198] 

that the exponent is insensitive to the corrugation in agreement with Monte-Carlo results [175]. The effect of 
the number of steps on the behavior of peff{z,t) for a given layer z, in the vicinity of the ordinary transition, 
is illustrated in Fig. 25. 

However, when the number of steps n is sufficiently smaller than the values of z, (3eff{z,n) is insensitive 
of the number of steps, when n increases more, (3eff{z,n) increases and once n exceeds the value of z, then 
(3eff{z,n) becomes independent of the number of steps, i.e after a certain degree of corrugation,/Je//(-2, ra) will 
be insensitive to the number of steps and then takes a saturation value, this value varies upon varying z. 

3.1.3 Wetting transition of a corrugated surface 

In the following, we denote by {Sk,k = 1, 2, ...) a configuration in which the top surface St and the first k — 1 
layers of the bulk have positive magnetizations; while the remaining M — k layers have negative ones. The 
ground state phase diagram is investigated exactly. In contrast to the case of perfect surface, in the case of 
short range substrate potential, the prewetting and layering transitions occurs at T = OK (Fig. 26). 

At finite temperature, phase diagrams obtained using Monte-Carlo simulations (Fig. 27) shows that the 
wetting temperature depends weakly on the surface corrugation degree but it depends strongly on lateral size 
of the system. Furthermore the behavior of the magnetization as a function of the external magnetic field is 
presented in Fig. 28a , for a fixed layer {z = 3) and several values of the position x in the x-direction. 

It is clear that, the increasing of the external bulk field, the top surface (x=18) exhibits a first order transition 
from negative to positive magnetization, accompanied by the transition of the bulk site located near the surface 
at z=3. This is due to both thermal fluctuations and the surface effect in the x-direction. The other bulk sites 
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Figure 26: Phase diagram of the ground state, in the [H / J. Hst/ J) plane, for a system with M — 10 layers, 
n = 2 steps and L — IQ spins per step, with the condition Hsb — —Hst [200]. 
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Figure 27: Layer transitions of the system, in the {H/ J, T/ J) plane, using Monte Carlo simulations for Hst / J ■ 
0.9 [200]. 
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Figure 28: Magnetisation profiles as a function of the reduced bulk field H/J for Hst/J = 0.9 and T/J — 3.5. 
The magnetizations m{x,z) are plotted for a fixed layer z = 3 and different positions x (a), and for a fixed 
position a; = and different layers z (b) [200] . 



keep their negative magnetizations. This is called "prelayer transition". Such behavior was obtained neither in 
the perfect surface nor at T = OK cases. However, by increasing H the so-called layer transition occurs when 
the magnetizations of the remaining bulk sites of the layer z,(z — 3) jump from negative values to positive ones. 
The effect of surfaces in the z-direction, is presented in Fig. 28b in which the local magnetization m{x, z) is 
presented as a function of H for a fixed value of x{x = 0) and several layers z. It is clear that the bulk field, at 
which the transition Sk ^ S'fc+i, (fc = 0, 1, M — 1) occurs, increases with increasing k. 

Moreover, using mean field theory, it is shown that the wetting temperature depends weakly on the corruga- 
tion degree n. Indeed, = 2.55 for n = 2, = 2.59 for n = 4 for a film with thickness N = 10. We note that 
the corrugate surface leads to the appearance of wetting in both longitudinal (z) and transverse (x) directions. 

3.2 Edge wetting in a three-dimensional system 

We consider a spin- 1/2 Ising model in a lattice formed with N layers, each one of them is constituted by two 
perpendicular perfect planes and contains 2(A'' — fc) -|- 1 spin chains which are infinite in the y-direction (Fig. 
29). The system is considered under the effects of an uniform surface magnetic field Hg applied on the planes 
{x, y,z = 1) ,{x = 1, y, z) of the layer k = 1 and an external magnetic field H applied on all spins of the system. 
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Figure 29: Geometry of the system formed with two surfaces {x,y,z = 1) and {x = l,y,z) with N spins in 
both the X and z— directions. The system is infinite in the direction. A uniform surface magnetic field 
is applied on the planes (z = l,x,y) and (x = l,y,z). An external magnetic field H is applied to the global 
system. 



However the Hamiltonian of the system is given by 

<i,j> i 

where, Si — ±1 are spin-1/2 Ising random variables, and = J are the exchange interactions assumed to be 
constant. Since the system is invariant by translation, in the following, the coordinate y will be cancelled. The 
total magnetic field Hi applied on a site 'i' of coordinates (x, z) is given by: 

{H + Hs for x=l,z = l,...,N 
H + Hs for z^l,x^l,...,N (35) 
H elsewhere 

H and Hs are, respectively, the external and surface magnetic fields both applied in the z-direction. 
The notation (Vlq) with p = 0,1, 2, N and q ^ 0,1, 3, 5, 2(A^ — p) — 1, will be used to denote that the 
first p layers and the q first spin chain of the layer p + 1 are in a magnetic state " up" ; while the remaining 
N — {p + 1) layers are in the state "down". In particular, the notation 1^ (resp. O^) will be used to denote 
a configuration with positive magnetization for all layer spins (resp. negative magnetization for all layer spins) 
of the system. Fig. 30 give a ketch of different possible configurations for a system of thickness = 4. 

Numerical results are given for two system sizes A^ = 4, using Monte Carlo simulations and A^ = 20 using 
mean field theory. However, the phase diagram obtained using the Monte Carlo (MC)simulations is represented 
in Fig. 31 for H^/J = 1.0. 

The behavior of the local magnetizations m(l,l), m(l,2), m(2,2) and m(2,3) are given in Figs. 32a and 
32b, as a function of the external magnetic field H . It is found that, when increasing the external field, the 
corners transits before the other intra-layering transitions. Such behavior was not obtained in the perfect 
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Figure 31: Phase diagram of the intra-layering and layer-by-laycr transitions, in the plane {H / J,T / J), using 
the Monte Carlo simulations for a system with N = A,ny = 100 and Hg/ J = 1 [40]. 



surface case[l-12,69,70]. However, the other layer-by-layer transitions occur when increasing the external field 
at specific temperature values. 

As one can expects, the external field values needed to make arising the intra-layering and inter-layering 
transitions [40] increase with increasing the order of the layer counted from the surface fc = 1. This is qualitatively 
in good agreement with previous works [1-12,69,70]. it is seen that the increasing system size effect is to decrease 
the intra and inter layering temperature values. The bulk layering transitions are shifted to higher external 
magnetic field values, for a fixed surface magnetic field value. 

The behavior of the surface, T£, and the bulk T| intra-layering temperature profiles as a function of the thickness 
N, for two values of Hs/J =1.0 and Hs/J = 0.9 is presented in Figs. 33. 

It is clear that T£ as well as T| decreases when increasing the system size. For a fixed system size N, these 
intra-layering temperatures decrease when increasing the surface magnetic field values. Moreover, reaches 
the wetting temperature for sufficiently large values of the system size. These results means that the edge leads 
to the appearance of the intralayering transitions which are absent in the case of perfect surfaces [1-12,45,69,70] 
and the continuous model (Milchev et al. [40]), [4 7-51]; On the other hand it leads to a decay of the wetting 
temperature compared with the perfect surface case [1-12,45,69,70]. This is in good agreement with the results 
obtained in Ref. [47-51]. 
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Figure 32: Magnetization profiles as a function of the reduced bulk magnetic field H/J for ./V = 4 and 
Hs/J = 1.0, of m(l, 1) and m(l, 2) at T/J = 0.75 (a); m(2, 2) and m(2, 3) at T/J = 1.0 (b), using the Monte 
Carlo method [40]. 




Figure 33: Surface layering temperature profiles as a function of the system size N for two surface magnetic 
field values Hs/J = 1.0 and Hs/J = 0.9, by using the mean field method [40]. 
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4 Quantum effects on wetting and layering transitions 



4.1 Uniform surface transverse field 

We consider N coupled ferromagnetic square layers in a transverse and longitudinal magnetic fields. The 
Hamiltonian of the system is given by 

H = -J2 J^oStS^i - J2(^St + HiSt) (36) 

<ij> i 

where, S", {a = x, z) arc the spin-1/2 Pauli matrices, is the transverse magnetic field in the x direction and 
Hi is the magnetic field in the z direction, defined by: 



Hi=l 



H + Hsi for i = l 

H for 1< i < TV (37) 

H + Hs2 for i = N 



The surface magnetic field Hs\ is applied on the first layer {i = 1), the last layer {i = N) is under a surface 
magnetic field Hs2 and H is the bulk magnetic field. 

For a simple cubic lattice, the magnetization and the free energy for each plane are given, respectively, by the 
following expressions: 

mj, = [4mp + r(mp+i + rup-i + H)] tanh(/JAp) (38) 



and 



with 



Fp = -i log(2cosh(/?Ap)) + ^nipiAnip + r{mp+i + nip-i) (39) 



Ap = ^ (4mp + r(mp+i + mp_i) + H)"^ + O^, (40) 

The state equations are solved numerically, it is found that there exist a wetting transverse field flw/J, above 
which a sequence of layering transitions occurs, which depends on the value of the temperature and the surface 
magnetic field, hsi- On the other hand, when T > T^, the layering transitions exist in the absence of the 
transverse magnetic field, this means that the temperature effect is sufficient to produce layering transitions. 
(T, O) Phase diagram is presented in Fig. 34. 

However, the wetting transverse magnetic field, flw, decreases with increasing temperature and/or decreasing 
surface magnetic field. 

4.2 Variable surface transverse field 

In this case we consider a surface variable transverse field, the transverse field of a layer z (distance from the 
substrate) is given by: 

= ft,/z" (41) 

with fig corresponding to the transverse magnetic field applied on the first layer z = I. Values of the parameter 
a will be discussed in the following, in particular, a = corresponds to a uniform transverse field applied over 
all the layers. The notation O^D^~p means that the top p layers are ordered while the remaining N — p are 
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Figure 34: Phase diagram in (T/J, 0/J) plane, for both cases: [Hgi/J = 0.8, Hs2/J = —Hgi/J) and {Hsi/J = 
0.9,Hs2/J = 0.), and N = 20 [79]. 



disordered. Phase diagrams for a = 1.0, is presented in Figs. 35. It is found that the layers with higher value of 
z need greater magnetic transverse field values to become disordered. 

However, for a fixed temperature, the disorder of each layer occurs when increasing the transverse magnetic 
field so that ^c{z = 1) < ilc{z = 2) < ... < flc{z = N). It is found that in absence of the transverse magnetic 
field, fig — 0, the critical temperature for all the layers, 1 < ^; < A'', is close to: Tj^^(bulk) = 6.0, for the mean 
field theory case, and T/^'"(bulk) = 5.073 for the finite cluster method. The corresponding thermal variation of 
the longitudinal and transverse magnetizations are computed using mean field theory and presented respectively 
in Fig. 36a and 36b. 

The behavior of the longitudinal magnetization critical exponent 7 at a fixed temperature T = 1.0 is also 
calculated as a function of the exponent a . However, 7 depends on the layer z, and decreases from 1, for small 
values of a, to 0.5 for higher values of a (Fig. 37) 

4.3 Random transverse field effect on order-disorder layering transitions 

Now we consider the case of a surface random transverse field governing by a two peaks distribution law namely: 

V{fli) = p6{fli - Qs) + (1 - p)S{fli) (42) 

whore p is the probability that Qi acts on a site i of the surface with the value fig. using finite cluster 
approximation within an expansion technique for spin-1/2 cluster identities, the longitudinal magnetizations 
m^, for a plane 'fc' (k = 1,...,N) are well calculated. However, On the other hand, it is shown that for a 
fixed a and very low temperatures, each layer k disorders at a critical probability Pc{k) corresponding to a 
critical surface transverse magnetic field O^(fc). Indeed, as it is shown in Fig. 38a, for a = 1.0, the film is 
ordered Sii p < Pc — 0.9949 even for very large values of the surface transverse magnetic field. While for 
Pc < P < Pc{k = 1) = 0.9950 (Fig. 38b), the first surface becomes disordered, increasing p, the second layer 
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Figure 35: The dependence of the critical temperature Tc on the surface transverse magnetic field, for a = 1.0 
(N=20 layers), using the mean field theory [253]. 




Figure 36: The dependence of the longitudinal layer magnetizations (a) as a function the transverse magnetic 
field for a = 1.0 and a fixed temperature: T = 1.0 ; and the thermal behavior of the transverse layer magneti- 
zations (6), for a = 1.0 and Og = 50.0. The number accompanying each curve denotes the value of the layer z 
[253]. 
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Figure 37: The dependence of the critical exponent 7 on the parameter a for T = 1.0. The number accompa- 
nying each curve denotes the value of the layer z, using the mean field theory [253]. 

disorders at -pcXk = 1) < p < Pc(fc = 2) = 0.9955 (Fig. 38c), and so on. Finally the film is totally disordered 
when the applied surface transverse magnetic field is acting on all the surface sites (i.e p = 1) (Fig. 38d). 

It arises that a sequence of order-disorder layering transitions occurs when decreasing probability from the 
value J) = 1.0, corresponding to a totally ordered film at very low temperatures. We note that the critical 
surface transverse field il^ depends on the value of a and the probability p. Indeed it is clear from Fig. 39 that 
rig increases when increasing a and/or decreasing p. 

On the other hand the profile of the surface critical transverse field (Fig. 40) and the surface critical 
probability (Fig. 41) shows that 17^ and increases with the depth of the film. 

The surface critical exponent of the longitudinal magnetization is also computed (Fig. 42) as a function of 
a for several values of p. It is found that ^{k) decreases when increasing a and/or decreasing -p. 

4.4 Random transverse field effect on wetting transitions 

Now we will consider the effect of random transverse field on the wetting transition of a spin-1/2 Ising fer- 
romagnetic film in both surface and external magnetic fields, described in Eqs.(36), (42). Using finite cluster 
approximation, numerical calculations are limited to a film with thickness. However, phase diagram in the space 
(T, f2,p) plotted in Fig. 43, shows the existence of a critical value Pcj of the probability p, above which the 
wetting transitions disappears. 

The value of depends on the value of Hg. However, Fig. 44 shows that Pc decreases when increasing 
Hg. indeed, decreases linearly for Hs < 0.85, and almost linearly for Hg > 0.85, with a discontinuity at 

= 0.85. 

A similar result was found by Harris [237], in the case of both bond and site dilutions, showing that the 
critical transverse field at zero temperature presents a discontinuity as the concentration passes through the 
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Figure 38: The dependence of the critical temperature Tc on the surface transverse magnetic field, for a = 1.0 
and different values of the probability: p{= 0.9945) < Pc{= 0.9949) (a), Pc < P < Pc{k = 1) = 0.9950 (b), 
Pcik = l)<p< Pc{k = 2) = 0.9955 (c) and p = 1.0 (d) [254]. 
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Figure 39: The dependence of the critical surface transverse magnetic field ^^(fc = 1) as a function of the 
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Figure 40: The dependence of the critical surface transverse magnetic field = 1) as a function of the 
number of layers for T = 0.05 and different values of the exponent a = 0.5, a = 1.0, a = 1.5 and a = 2.0 [247]. 
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Figure 42: The dependence of the surface critical exponent 7(fc = 1) on the parameter a for a fixed temperature 
T = 1.0 and different values of the probability p [247]. 
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Figure 44: The critical probability pc dependence on the surface field Hgi/J for T/J — 0.02 [247]. 
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critical percolation concentration. 



4.5 Layering sublimation transitions of the Blume-Emery-GrifRths in a transverse 
field 

In this part, we consider a spin-1 BEG model on a simple cubic lattice, governed by the Hamiltonian 

H = -JY^ S^.S,. -KJ2 SlS% -Aj2sl-nY,S,, (43) 
{id) (id) i i 

where J and K represent respectively the bilinear and bi-quadratic couplings, A is the chemical potential in the 

lattice-gas interpretation [246,249], fl is the transverse field, indicates summation over nearest-neighbor 

sites of the simple cubic lattice, and Six, Using the mean field theory, the free energy, the longitudinal 

magnetization and quadrupolar moments are given respectively by: 

TV 3 

f = Y. -KbTLoqIY, exp{-(3Xk{p))] + -m,(p)(TO,(p - 1) + 4m,(p) + m,(p + 1)) (44) 

p=l k=l 

+fQz{p){Qz{p - 1) + Mz{p) +qz{p+ 1)) 



ELiK(p)-/?g(p))exp(-/3Ai.(p)) 
Ek=iexp(-/3Ak(p)) 



^^(p) ^ ELiK(p)+/3^(p))exp(-/?A,(p)) ^^g^ 
Ek=i exp(-/3Ak(p)) 



where 



with 



\k{p) = \{D{p) + y^cos(^fe(p))) (47) 

a iv) = \il{\k{p)-D{p)+A{p))\ 

2VV02((Afc(p) - D{p)f + A^p)) + ((Afc(p) - D{p)f - A^{j>)f 

Xkjp) - D{p) - A{p) 

Mp) = ^arccos{-27q(p)/2Q{p)) + 2/3(fc - l)7r (50) 

Q{P) = ^y2VMTMpy+2VM (51) 

r{p) = -{A^{p) + n')-D^p)/3 (52) 

q{p) = -D{p){2A\p) - {2/9)D\p) - n^)/3 (53) 

A{p) = -J{m,{p - 1) + Am,{p) +m,{p+ 1)) (54) 

D{p) = - A - K{qz{p - 1) + 4(7, (p) +qz{p+ 1)) (55) 
With the free boundary conditions, mz{N + 1) = TOz(O) = qz{N + 1) = (?z(0) — 0. 

Equations (8) and (9) are the order parameters of the system which enable us to characterize transitions. Their 
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Figure 45: Phase diagram in the (O, A) plane for T = [248]. 



calculations have been performed numerically. However, in the case of several distinct solutions. The one which 
have the smaller free energy (19) is accepted as the physical solution. However, numerical calculations are 
done for N = 64 and K/J = 3 because at the latter value and in the absence of the transverse field, phase 
diagram of the spin-1 BEG model in temperature and chemical potential plane is equivalent to the solid-gas- 
liquid phase diagram [246]. The solid, liquid and vapor phases are characterized by a high density and non-zero 
magnetization, a high density and zero magnetization and a low density and zero magnetization respectively 
[246]. We note that the bulk solid- vapor coexistence take place at a chemical potential Acoe = — 12 + 4:/N in 
the case T = and O = 0. Hence for N=64, A^oe = -11.9375. Hence, we denote here after by y2m5iv-2m ^ 
m = 0, 1, N/2 a configuration of the system with m top layers and m bottom layers in the vapor state, while 
the remaining layers of the system are in the solid state; and by L'^ a configuration with N layers in the liquid 
state. and are respectively the bulk vapor and solid states respectively. 

Phase diagram at T=Ok (Fig. 45) shows that a sequence of layering sublimation transition (y2m^iv-2m ^ 
Y2m+2gN-2m-2^ appears at diflFerent values of the transverse field Clm for a fixed value of chemical potential 
A) — 11.9375 (i.e in the bulk solid region); and then the bulk melting transition (5^ ^ L^) occurs at a critical 
bulk transverse field fisL- 

While for A(Acoe and 0(0c there is also the bulk vapor-liquid transition {V^ L^) which take place at a 
critical value of a transverse field O-vl, while the solid- liquid-gas bulk coexistence occurs at the triple transverse 
field Q.T and A = Acoe- Such result is qualitatively similar to the one obtained by Gelb [246] under the effect of 
temperature. In Fig. 46, the bulk density exhibits a discontinuity at the transition low density (Vapor phase) 
to high density (liquid phase) which occurs at Q.vl = 5.51 for A = —12.002. 

The behavior of density and magnetization, as a function of the transverse field f2, are given, for A)Acoe, in 
Fig. 47a and 47b respectively. It is clear that the transition of a layer p (from solid state to quasi-vapor state) 
is characterized by the jump of both density and magnetization. 

At T > 0, the (T, O) phase diagrams presented in Fig. 48 show the existence of a sequence of layering 
sublimation transition for each value of temperature and fixed value of chemical potential, and a solid-liquid 
bulk transition occurs at a transverse field Q.sl which depends on the value of chemical potential. 

However Q,sl increases with increasing chemical potential. For example for T = 1 and A = —11, O-sl = 6.01 
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Figure 46: The dependence of the bulk density of the transverse field for A = 12.002 and T = 0. The vapor 
linnid transition is found at Qi/r = 5.51 [2481. 
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Figure 47: (a) Profile of density of the top three and last layers as a function transverse field for A = 11 
and T = 0. (b): Profile of magnetization of the first three and last layers as a function of transverse field for 
A = -11 and T = 0K [248]. 
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(Fig. 48a) while it takes the value 5.30 for T=l and A = —11.9 (Fig. 48b). On the other hand the transition line 
of a layer p terminates by an end point of coordinate (Tc^, ficp)- these endpoints converge towards a roughening 
temperature Tr and roughening transverse field CIr above which all interfaces are rough on a microscopic scale, 
Hence Tr — 3.07 and Qr = 5.18 in the case of Fig. 48a. Besides, in this figure the critical end point of the 
first layer is greater than the second layer end point while in Fig. 48b the critical end point of the first layer 
is smaller than the second layer one. This means that, depending on the value of the chemical potential, there 
exist a critical temperature Tg and critical transverse field Clg at which the surface transition occurs while the 
other layers transitions cannot take place for any value of transverse field. 
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